電磁場のディラック方程式（maxima ソース）
電磁場のディラック方程式の計算に使った maxima (wxMaxima) のソースです。コメント等ありませんが，だいたい見たらわかると思います。gt 〜 gz がγ行列です。
電流項になにか計算ミスがあるので L と L2 のふたつ行列を定義しましたが，本当は同じはずです。
中嶋さんのご指摘により解決しました。(2019/4/11)
code:maxima
* wxMaxima batch file version 1 DO NOT EDIT BY HAND! */
/* Created with wxMaxima version 13.04.2 */
/* wxMaxima: input start */
g0:;
/* wxMaxima: input end */
/* wxMaxima: input start */
gt: matrix(
0,0,1,0,
0,0,0,1,
1,0,0,0,
0,1,0,0
);
/* wxMaxima: input end */
/* wxMaxima: input start */
gx: matrix(
0,0,0,1,
0,0,1,0,
0,-1,0,0,
-1,0,0,0
);
/* wxMaxima: input end */
/* wxMaxima: input start */
gy: matrix(
0,0,0,-%i,
0,0,%i,0,
0,%i,0,0,
-%i,0,0,0
);
/* wxMaxima: input end */
/* wxMaxima: input start */
gz: matrix(
0,0,1,0,
0,0,0,-1,
-1,0,0,0,
0,1,0,0
);
/* wxMaxima: input end */
/* wxMaxima: input start */
D:(Dt*gt-Dx*gx-Dy*gy-Dz*gz);
/* wxMaxima: input end */
/* wxMaxima: input start */
F: (Ex* gt.gx+ Ey* gt.gy + Ez*gt.gz) + (Bz*gx.gy+By*gz.gx+Bx*gy.gz) ;
/* wxMaxima: input end */
/* wxMaxima: input start */
J:(Jt*gt-Jx*gx-Jy*gy-Jz*gz);
/* wxMaxima: input end */
/* wxMaxima: input start */
K: matrix(
1,1,1,1,
1,-1,1,-1,
1,1,-1,-1,
1,-1,-1,1
);
/* wxMaxima: input end */
/* wxMaxima: input start */
L:matrix(1,1,1,1)/4;
/* wxMaxima: input end */
/* wxMaxima: input start */
realpart(expand(K.D.F.L));
/* wxMaxima: input end */
/* wxMaxima: input start */
imagpart(expand(K.D.F.L));
/* wxMaxima: input end */
/* wxMaxima: input start */
realpart(expand(K.J.L));
/* wxMaxima: input end */
/* wxMaxima: input start */
imagpart(expand(K.J.L));
/* wxMaxima: input end */
/* Maxima can't load/batch files which end with a comment! */
"Created with wxMaxima"$
＝＝＝＝＝＝＝＝＝＝＝＝＝＝
#クリフォード代数
#電磁場のディラック方程式
このサイトについて